Stability of Quantum Critical Points in the Presence of Competing Orders 
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We investigate the stability of Quantum Critical Points (QCPs) in the presence of two competing 
phases. These phases near QCPs are assumed to be either classical or quantum and assumed to 
repulsively interact via square-square interactions. We find that for any dynamical exponents and 
for any dimensionality strong enough interaction renders QCPs unstable, and drives transitions to 
become first order. We propose that this instability and the onset of first-order transitions lead to 
spatially inhomogeneous states in practical materials near putative QCPs. Our analysis also leads 
us to suggest that there is a breakdown of Conformal Field Theory (CFT) scaling in the Anti de 
Sitter models, and in fact these models contain first-order transitions in the strong coupling limit. 



INTRODUCTION 



^ Quantum criticality is an important concept that has dominated the landscape of modern condensed matter physics 

for the last decade The idea behind quantum criticality is simple and powerful. Imagine competing interactions 
that typically drive the transitions between different phases. Logically one has to allow for the possibility that the 
relative strength of these competing interactions is tunable as a function of the external control parameters such 
as pressure, magnetic field, doping: we deliberately omit temperature as a control parameter since quantum phase 
transitions (QPTs) will occur at T=0. The simplest route to arrive at a QPT is to consider a line of finite temperature 
. phase transition as a function of some control parameter, such as pressure P, magnetic field B or doping x. At T = 
this line will indicate a critical value of the control parameter. This specific value of the control parameter, where one 
expects a precise balance between tendency to different phases or states, is called a quantum critical point (QCP). 
Near this point, competing interactions nearly compensate each other. It is often asserted that it is the physics 
of frustration and competition, which leads to the finite temperature transition, that also controls and enables the 
interesting properties of materials that are brought to the T = QCP. 

Much of the attention on quantum criticality has been focused on the finite temperature scaling properties 1-3]. 
Temperature is the only relevant scale in the quantum critical region above the QCP, bounded by the crossover line 
OO ■ T* ~ \r\ uz . The parameter r measures the distance to the QCP, v is the correlation-length exponent in £ ~ r" 
t-H ' and z is the dynamical exponent in £ T ~ £ z . With the correlation length £ and correlation time £ T much larger 
, than any other scale of the system, power law behavior is expected for many physical observables, e. g. the specific 
heat, magnetic susceptibility, and most notably resistivity. Clear deviations from the Fermi liquid predictions are 
experimentally detected, and these phases are commonly termed no n- Fermi liquids. In many systems, the anomalous 
finite temperature scaling properties are asserted to result from the underlying zero temperature QCPs. 
^ , In this paper, we would like to emphasize another aspect of quantum criticality, namely that it serves as a driving 
force for new exotic phenomena at extremely low temperatures and in extremely clean systems. One possibility is the 
appearance of new phases around the QCPs. It has been found in numerous experiments as one lowers temperature, 
seemingly inevitably in all the systems available, new phases appear near the QCP. Most commonly observed to date 
is the superconducting phase. The phenomenon of a superconducting dome enclosing the region near the QCP is quite 
general (see Fig. 1). It has been identified in many heavy fermion systems 4-6], plausibly also in cuprates [7j, even 

1 6|. O ther examples include the nematic 
20], the origin of which is still under 



possibly in pnictides [l-Qjl , and probably in organic charge-transfer salts 1 




phase around the metamagnctic QCP in the bilayer ruthenate S^R^Oy 

intense debate fill-Eol]. The emerging quantum paraelectric - ferroelectric phase diagram is also very reminiscent 
27], as is the disproportionation-supcrconducting phase in doped bismuth oxide superconductors [28l - l33j . 
It has also been discovered recently that, as samples are becoming cleaner, on the approach to QCP we encounter first 
order transitions, and the new phases near the QCP are usually inhomogeneous and exhibit finite wavevector orderings 



(see 3_4j-|36| and references therein). For example, the heavy fermion compound CeRhlns orders antiferromagnetically 
at low temperature and ambient pressure. As pressure increases, the Neel temperature decreases and at some pressure 
the antiferromagnetic phase is replaced by a superconducting phase through a first-order phase transition. There are 
also evidences for a competitive coexistence of the two phases within the antiferromagnetic phase, as in some organic 



2 




FIG. 1: Illustration of the competing phases and superconducting dome. Here for concreteness, we consider the ordered phase 
to be an antiferromagnetic phase, x is the tuning parameter. It can be pressure, magnetic field or doping. The superconducting 
temperature usually has the highest value right above the QCP. 



charge-transfer superconductor precursor antiferromagnetic phases. Such coexistence was also observed in Rh-doped 
Celrlns . The heavy fermion superconductor CeCoIns has the unusual property that when a magnetic field is applied 
to suppress superconductivity, the superconducting phase transition becomes first-order below To ~ Q.IK . For the 
superconducting ferromagnet UGe2, where superconductivity exists within the ferromagnetic state, the two magnetic 
transitions (ferromagnetic to paramagnetic and large-moment ferromagnetic to small-moment ferromagnetic) are both 
first order I37H 391 . Other examples of continuous phase transitions turning first-order at low temperatures include 



CeRh 2 Si 2 |4(J,|4l|, Celn 3 [42j, URhGe [43J, ZrZn 2 |44j and MnSi [45[. The prevailing point of view seems to be that 
this happens only in a few cases and these are considered exceptions. Yet we are facing a rapidly growing list of these 
"exceptions", and we take the view here that they rather represent a general property of QCPs. 

The point is that, on approach to the QCP, an interaction that was deemed irrelevant initially, takes over and 
dominates. For example it has been proposed recently that the superconducting instability, which is marginal in the 
usual Fermi liquids, becomes relevant near the QCP and leads to a high transition temperature (46|. Actually these 
instabilities are numerous and can vary, depending on the system at hand. However there seems to be a unifying theme 
of those instabilities. We suggest that QCPs are unstable precisely for the reasons we are interested in these points: 
extreme softness and extreme susceptibility of the system in the vicinity of QCPs. We regard the recently discovered 
first order transitions as indicators of a more fundamental and thus powerful phsyics. We are often prevented from 
reaching quantum criticality, and often the destruction is relatively trivial and certainly not as appealing and elegant 
as quantum criticality. We can draw an analogy from gravitational physics, where the naked singularities are believed 
to be prevented from h app ening due to many kinds of relevant instabilities. This is generally known as the "cosmic 
censorship conjecture" (47]] ■ The recently proposed AdS/CFT correspondence 48- 5(|, which maps a non-gravitational 
field theory to a higher dimensional gravitational theory, adds more to this story. Here researchers have begun to 
realize that the Reissner-Nordstrom black holes in AdS space, which should have a macroscopic entropy at zero 
temperature, are unstable to the spontaneous creation of particle-antiparticle pairs, and tend to collapse to a state 
with lower entropy 51 - 53| . 

There have appeared in the literature scattered examples of first-order quantum phase transitions at the supposed-to- 
be continuous QCPs [2, 54 -6(| , however it appears that the universality of this phenomenon is not widely appreciated. 
This universality is the main motivation for our paper. We will systematically study the different possibilities for 
converting a continuous QPT to first order. 

The first striking example how fluctuations of one of the order parameters can qualitatively change the nature 



of the transition comes from the Coleman- Weinberg model [6l|, where they showed how gauge fluctuations of the 
charged field introduce a first order transition. In this work it was shown that in dimension d = 3, for any weak 
coupling strength, one develops a logarithmic singularity, and therefore the effective field theory has a first-order phase 
transition. Subsequently, this result was extended to include classical gauge field fluctuations by Halperin, Lubensky 
and Ma [r32j |. where a cubic correction to the free energy was found. Nontrivial gradient terms can also induce an 
inhomogeneous phase and/or glassy behavior (r33l |. 

A prototypical example for the competing phases and superconducting dome is shown schematically in Fig.l. 
Below, we apply the renormalization group (RG) and scaling analysis to infer the stability if the QCP as a result 
of competition. We find in our analysis that the QCP is indeed unstable towards a first order transition as a result 
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of competition. Obviously details of the collapse of a QCP and the resulting phase diagram depend on details of 
the nature of the fluctuating held and details of the interactions. We find that the most relevant parameters that 
enter into criterion for stability of a QCP are the strength of interactions between competing phases: we take this 
interaction to be repulsive between squares of the competing order parameters. When the two order parameters break 
different symmetries, the coupling will be between the squares of them. Another important factor that controls the 
phase diagram is the dynamical exponents z of the fields. The nature of the competition also depends on the classical 
or quantum character of the fields. Here by classical we do not necessarily mean a finite temperature phase transition, 
but rather that the typical energy scale is above the ultraviolet cutoff, and the finite frequency modes of the order 
parameters can be ignored, so that a simple description in terms of free energy is enough to capture the physics. We 
analyzed three possibilities for the competing orders: 

i) classical + classical. Here we found that interactions generally reduce the region of coexistence, and when 
interaction strength exceeds some critical value, the second-order phase transitions become first order. 

ii) classical + quantum. Here the quantum field is integrated out, giving rise to a correction to the effective potential 
of the classical order parameter. For a massive fluctuating field with d + z ^ 6, or a massless one with d + z ^ 4, the 
second-order quantum phase transition becomes first order. 

iii) quantum + quantum. Here RG analysis was employed, and we found that in the high dimensional parameter 
space, there are generally regions with runaway flow, indicating a first-order quantum phase transition. 

It has been proposed recently that alternative route to the breakdown of quantum criticality is through the basic 
collapse of Landau- Wilson paradigm of conventional order parameters and formation of the deconfined quantum 
critical phases ([ill, EH). This is a possibility that has been discussed for specific models and requires a different 
approach than the one taken here. We are not addressing this possibility. 

The plan of the paper is as follows. In section II, we consider coupling two classical order parameter fields together. 
Both fields are characterized by their free energies and Landau mean field theory will be used. In section III, 
we consider coupling a classical order parameter to a quantum mechanical one, which can have different dynamical 
exponents. The classical field is described by its free energy and the quantum held by its action; the latter is integrated 
out to produce a correction to the effective potential for the former. In section IV, we consider coupling two quantum 
mechanical fields together. With both fields described by their actions, we use RG equations to examine the stability 
conditions. In particular, we study in detail the case where the two coupled order parameters have different dynamical 
exponents, which, to our knowledge, has not been considered previously. In the conclusion section, we summarize our 
findings. 

TWO COMPETING CLASSICAL FIELDS 

We consider in this section two competing classical helds. Examples are the superconducting order and antiferro- 
magnetic order in CeRhlns and Rh-doped Celrlns, and the superconducting order and ferromagnetic order near the 
large-moment to small-moment transition in UGc2. We will follow the standard textbook approach, and this case is 
presented as a template for the more complex problems studied later on. 

We first study the problem at zero temperature. For simplicity, both of them are assumed to be real scalars. The 
free energy of the system consists of three parts, the two free parts i 7 ^, Fm and the interacting part F{ nt : 

F =F^ + Fm + F int ; 

n^(v^) 2 -aV 2 + §V> 4 ; 

a 8 (1) 

F M = P -f{VMf - a M M^ +PfM i - 

F int ^ 2 M 2 . 

Here, by changing a, olm, the system is tuned through the phase transition points. When the two fields are decoupled, 
with 7 = 0, there will be two separated second-order phase transitions. Assume the corresponding values of the 
tuning parameter x at these two transition points are x\ and X2, we can parameterize a, cum as a = a{x — x\) and 
ctM = o-m(x2 — %), where a, om are constants. 

We would like to know the ground state of the system. Following the standard procedure, we first find the 
homogeneous field configurations satisfying ^ = — 0, and then compare the corresponding free energy. It is 

easy to see that the above equations have four solutions, with (|^|, \M\) = (0, 0), (0, \Jcim / Pm), (y/a//3, 0), M*), 
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FIG. 2: (Color online) Illustration of the mean field phase diagram for two competing orders. Here for concreteness we consider 
antiferromagnetic and superconducting orders. The two orders coexist in the yellow region, whose area shrinks as the coupling 
increases from left to right. The left figure has 7 = 0, the central one has < 7 < \fppM, and the right one has 7 > ^//3/3a/- 
When 7 exceeds the critical value ^/PPm, the two second-order phase transition lines merge and become first order (the thick 
vertical line). 
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and the rescaled parameters are 7' = ^/aa Ml (3' = P/o?,p M — (3 M /a 2 M . When 7 = 0, the fourth solution reduces to 
(^*,M«) = (y/a/P, y/ani / Pm), with the two orders coexisting but decoupled. We are interested in the case where 
the two orders are competing, thus a relatively large positive 7. 

For x\ < x < X2, we have a > 0, ccm > 0. The necessary condition for the existence of the fourth solution 
is 7' > p, j3' M , \J P'P' M or 7' < p, p M , \J PP M - In this case, the configuration (0,0) has the highest free energy 
F[0, 0] = 0. For the configuration (i/j*, M*) with coexisting orders to have lower free energy than the two configurations 
with single order, one needs to have 7' < y/PP M , which reflects the simple fact that when the competition between 
the two orders is too large, their coexistence is not favored. Thus the condition for the configuration (ip*, M*) to be 
the ground state of the system is 7' < p and 7' < p M . If 7' > min{/3', P' M }, one of the fields has to vanish. 

Next we observe that, for x near Xi, p M remains finite, a ~ (x — xi), and 7' diverges as l/(x — xi), while p diverges 
as l/(x — x\) 2 . So the lowest energy configuration is ip = 0, \M\ = \J olm / Pm ■ Similarly, near #2, the ground state is 
(y/ a/ P 0). The region with coexisting orders shrinks to 



-fQM + pMCl 



- /3 M axi jaxi + BaMXi 

< x < 



7a + /3ttM 



(3) 



For 7 < vppjvf, this region has finite width. In this region, (0, 0) is the global maximum of the free energy, 
(0, \/ctM / Pm), (y/a/P,0) are saddle points, and (?/>», Af*) is the global minimum. The phase with coexisting or- 
der is sandwiched between the two singly ordered phases, and the two phase transitions are both second-order. The 
shift in spin-density wave ordering and Ising-nematic ordering due to a nearby competing superconducting order has 
been studied recently by Moon and Sachdev (66l . \6l\ , where they found that the fermionic degrees of freedom can play 
important roles. 

For 7 > \/PpMi this intermediate region with coexisting orders vanishes, and the two singly ordered phases are 
separated by a first-order quantum phase transition. The location of the phase transition point is determined by 
equating the two free energies at this point, 



F 



la(x c ) 



,0 



F 



0, 



P. 



M 



(4) 



which gives x c = (x 2 + Axi)/(l + A), with A = (aj o,m)\J Pm I P- The slope of the free energy changes discontinuously 
across the phase transition point, with a jump 



SF^ = 



dF 

fix 



dF\ 
dx) ; 



aau 



VPP 



(x 2 - Xi). 



(5) 
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The size of a first-order thermal phase transition can be characterized by the ratio of latent heat to the jump in 
specific heat in a reference second-order phase transition [62j . A similar quantity can be defined for a quantum phase 
transition, where the role of temperature is now played by the tuning parameter x. We choose as our reference point 
7 = 0, where the two order parameters are decoupled. For x < X\, one has d 2 F/dx 2 = — a 2 M / Pm\ f° r x > x 2> one has 
d 2 F/dx 2 — —a 2 / (3; and d 2 F/dx 2 = —cl 2 m /I3m — a? / (3 for xi < x < x%. We take the average of the absolute value of 
the two jumps to obtain 

6FW = l -{a 2 M /p M + a 2 /P). (6) 
So the size of this first-order quantum phase transition is 

6x= wm = WJ -(x 2 -x 1 ), (7) 

with j3 = f3/a 2 and /3m — /3m/c^m- It is of order X2 — x\, when /3 and (3m are not hugely different. 

The above consideration can be generalized to finite temperature, by including the temperature dependence of all 
the parameters. Specially, there exists some temperature T*, where x\{T*) = 22 (T*). In this way we obtain phase 
diagrams similar to those observed in experiments (see Fig. 1). 



EFFECTS OF QUANTUM FLUCTUATIONS 

In this section, we consider coupling an order parameter ip to another field <fi, which is fluctuating quantum me- 
chanically. The original field if) is still treated classically, meaning any finite frequency modes are ignored. For the 
quantum fields, in the spirit of Hertz-Millis-Moriya [68|-|70j], we assume that the fermionic degrees of freedom can be 
integrated out, and we will only deal with the bosonic order parameters. This model may, for example, explain the 
first-order ferromagnetic to paramagnetic transition in UGe2 , where the quantum fluctuations of the superconducting 
order parameter are coupled with the ferromagnetic order parameter, which can be regarded as classical near the 
superconducting transition point. 

We will integrate out the quantum field to obtain the effective free energy of a classical field. The partition function 
has the form 

Z[il>(r)] = J Vcf>(r, t) exp (-^ - S 4 - . (8) 

The free energy is of the same form as in the previous section with J-^ = j d d rF^. Thus, in the absence of coupling 
to other fields, the system goes through a second-order quantum phase transition as one tunes the control parameter 
x across its critical value. We consider a simple coupling 

Sw=g J d d rdT^j 2 <j) 2 . (9) 

The action of the cf> field depends on its dynamical exponent z. We notice that such classical + quantum formalism 
has been used to investigate the competing orders in cuprates in [7l| . 
The saddle point equation for ip reads 



6lnZ[ip(r) 
<5V>(r) 

which gives 



0, (10) 



-a + /^ 2 (r) - |V 2 + g(4> 2 (r))] ^(r) = 0. (11) 
Here we have defined the expectation value, 

(0 2 (r)) = ~y P^(r',T')| ^ 2 (r,r)exp(-5 -^). (12) 
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It can also be written in terms of the different frequency modes, 

(0 2 (r)> =T^(0(r, Wn )0(r,- W „)) =t£ f V^r 1 , u s )4>{r, tj n )cf>(r, -u n ) exp (S^ - S^) . (13) 

The quadratic term in S<f> is of the form 

S f =E/ dV / rfVXr'^^Xo 1 ^^",^)^",-^), (14) 
or more conveniently, in terms of momentum and frequency, 

s f } =£ / ^^MOx^CM.M-k,-*.)- (15) 

So in the presence of translational symmetry, we find 

2 _ r d d k 1 

This leads to the 1-loop correction to the effective potential for -0, determined by 



(16) 



Sip 



2g(<f)il>. (17) 



So far we have been general in this analysis. Further analysis requires us to make more specific assumptions about 
the dimensionality and dynamical exponents. 

Fluctuations with d — 3, z — 1 

When the <f> field has dynamical exponent z = 1, its propagator is of the form 

Xo(k,u Jn ) = - rT ^ Tr ,. (18) 

A special case is a gauge boson, which has zero bare mass, and thus £ —> oo. This problem has been studied in detail 
by Halperin, Lubensky and Ma [(32| for a classical phase transition (see also (z3), and by Coleman and Weinberg 
[6l| for relativistic quantum field theory. Other examples are critical fluctuations associated with spin-density wave 
transitions and superconducting transitions in clean systems. We also note that Continentino and collaborators have 
used the method of effective potential to investigate some special examples of the fluctuation-induced first order 
quantum phase transition [2|, l55l - l58l | . 

Let us consider T — 0, for which the summation T^2 U can be replaced by the integral / doj/(2ir). We then get 
for the one-loop correction to the effective potential 

^"M-w/g^ ..., i (19) 



dip " T J 2w J (27r) d lu 2 + k 2 + + gf 2 ' 

Carrying out the frequency integral, we obtain for d = 3, 



nffW _ 9t f\ k p (20) 



H 2tt 2 J q ^k 2 + £- 2 + g*P 2 ' 

where an ultraviolet cutoff is imposed. Integrating out momentum gives 
5V$m _ g^P 



Sip 47T 2 



a VA^TF^ - (r 2 + g^) m WA^ + # 2 



(21) 
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FIG. 3: (Color online) Schematic illustration of the fluctuation-induced first-order phase transition. Here, for concreteness, we 
consider ferromagnetic and superconducting orders. The ferromagnetic order is regarded as classical, while the superconducting 
one as quantum mechanical. At low temperatures, the second-order ferromagnetic to paramagnetic phase transition becomes 
first order (the thick vertical line), due to fluctuations of the superconducting order parameter. 



which can be simplified as 

Sip 47T 2 

Combined with the bare part, 



A 2 + i(r 2 +^ 2 )-(r 2 +^ 2 )in 



2A 



yr 2 +^ 2 



we get the effective potential to one-loop order, 



167P 



-(r 2 +sv> 2 ) 2 m 



2A 



V^ 2 + 



(22) 



(23) 



(24) 



with the quadratic and quartic terms renormalized by a — a — <?(4A 2 + £~ 2 )/(327r 2 ) and /3 = (3 + 3g/(327r 2 ). When 
<f> field is critical with £ — > oo, the third term is of the well-known Coleman- Weinberg form ip A ln(2 A / \J gip 2 ) , which 
drives the second-order quantum phase transition to first order. 

For £ large but finite, we can expand the third term as a power series in £,~ 2 /(gip 2 ), and the effective potential is 
of the form 



1 



' (2r 2 # 2 + .gV)in 2A 



gi> 2 



(25) 



In addition to the Coleman- Weinberg term, there is another term of the form ip 2 \nip, and again we have also a 
first-order phase transition. 

To study the generic case where the cf> field is massive, we rescale the ip field and cutoff, defining 



u 2 = 



g^ 2 



A = 



2A 

F 1 ' 



The rescaled effective potential takes the form 

V cS (u) 

which can be further simplified as 



Au 2 + ^Bu 4 -(l + u 2 ) 2 ln 



A 



(26) 



(27) 



Vcs(u) = -Au 2 + -Bu 4 + (1 + u 2 ) 2 ln(l + u 2 ). 



(28) 
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FIG. 4: (Color online) The effective potential as a function of the rescaled field u for various parameters in the case d = 3, z = 1. 
Here V c s(u) = -Au 2 + ±Bu 4 + (1 + u 2 ) 2 ln(l + u 2 ), with B = -5, and A = -0.25, -0.116, from top to bottom. 



The above potential is plotted in Fig. 4. We notice that with large enough cutoff A, one generally has B = B — In A 
large and negative. For A < 1, u = is a local minimum. There are also another two local minima with u = y a 
positive solution of equation 

2(l + j/)ln(l+2/) + (l + B)j/ + l- A = 0. (29) 
So we generally have a first-order quantum phase transition in this case (see Fig. 3 for a schematic picture). 



Fluctuations with d = 3, z = 2 



With dynamical exponent 2 = 2, the propagator of field is 



Xo(k,w n ) 



\uj n \T + k 2 + t 



(30) 



Examples are charge-density-wave and antiferromagnetic fluctuations. In the presence of dissipation, superconduct- 
ing transitions also have dynamical exponent z = 2. So the one- loop correction to the effective potential at zero 
temperature becomes 



Sip 



, duj f d d k 
2# / ^ 



2tt J (2ir) d |w|t + fc 2 + £~ 2 + gip 2 



(31) 



The momentum integral is cutoff at |k| = A, and correspondingly the frequency integral is cutoff at \oj\tq = A 2 . First, 
we integrate out frequency to obtain 



-d) 

eg 

Sip 



7T J T Jo 



dkk 2 In f 



A 2 



k 2 + £~ 2 + gip 2 



(32) 



and then integrate out momentum, with the final result 



8V$[hli\ _ giP 
Sip 3tt 3 tq 



A 3 In 



£~ 2 + gip 2 + 2A : 



£- 2 + gip 2 + A 2 



+ 2A 3 + 2(C 2 + gijj 2 ) 3/2 arctan 



A 



-2(C 2 + gip 2 + A 2 ) 3/2 arctan 



vt- 2 + vV 

A 



v/£" 2 + gip 2 + A 2 



(33) 



Up to order A , this is 



Sip 



37T 3 To 



A 3 (2 + ln2 - |) + ^A(r 2 + # 2 ) + *(C 2 + 9i> 2 ) 



2\3/2 



(34) 




FIG. 5: (Color online) The effective potential as a function of the rescaled field u for (a) d — 3, z = 2, where we have plotted 
V c a(u) = -Au 2 + \Bu A + (1 + u 2 f /2 - 1, with B = -8, and A = -3, -2.597, -2.2 from top to bottom; (b) d = 3, z = 3, where 
V c ft(u) = -Au 2 + \Bu A + (l + u 2 f\n{l + u 2 ), with B = -10, and A = 0.1,0.208,0.3 from top to bottom; (c) d = 1, z = 2, where 
we have plotted V cff (u) = -Au 2 + \Bu 4 - (1 + w 2 ) 3/2 + 1, with B = 0.1, and A = -5.3, -5.1413, -5 from top to bottom; (d) 
d=l,z = l, where we have plotted V eB (u) = -Au 2 + ±Bw 4 - (1 + u 2 ) ln(l + u 2 ), with B = 0.3, and A = -1.45, -1.412, -1.39 
from top to bottom. All these plots are of similar shape. However, we notice that the scales are quite different. 



The first two terms just renormalize the bare a and fi. When the <fi field is critical, £ — > oo, the third term becomes 
of order ip 5 , and is thus irrelevant. When £ is large but not infinite, we get the effective potential 



2 lbn^To lbn^To 



In addition to the ip 5 term there is another term of order ip 3 , which may drive the second-order quantum phase 
transition to first order. 

Let us consider a massive <p field. Carrying out the same rescaling as we made for 2 = 1, we get the rescaled effective 
potential of the form 

V cS (u) = -Au 2 + X -Bu A + (1 + u 2 fl 2 . (36) 
For large negative B, we obtain a first-order quantum phase transition (see Fig. 5(a)). 



Fluctuations with d = 3, z = 3 



When the <f> field has dynamical exponent z — 3, e.g. for ferromagnetic fluctuations, its propagator is 

Xo(k,^„)= |*„| ,! 2 ,,_ 2 - ( 3? ) 

Thus the one-loop correction to the effective potential at T = is determined from 

SV^lip] „ , f dcu [ d d k 1 



^ 77" ,o.w ... (38) 



<5V 7 2n J (2ir) d 7 M + jfc2 + £-2 + ^ 
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with a momentum cutoff at |k| = A, and a frequency cutoff at 7|w| = A 3 . The frequency integral gives 



f 

Jo 



Stp 47r 4 7 
and the momentum integral further leads to the result 



dkk 3 In 



A 3 



fc 3 + fc(^- 2 + gV) 



aip 2 + \h a + ^^(r 2 + g^ 2 f in(r 2 + # 2 )- 



(39) 



(40) 



When <j> is critical, £ — > oo, the third term is of the form ip & ln^, which is irrelevant. For finite £, there is also a term 
of the form ip 4 ln?/>, which will drive the second-order quantum phase transition to first order. 
For general £, the rescalcd effective potential reads 



1 



Au 2 + -Bu 4 + (1 + u 2 f ln(l + u 2 ), 



V cS (u) 

which can lead to a first-order quantum phase transition (see fig. 5(b)). 



(41) 



Fluctuations with d = 3, z = 4 

For a dirty metallic ferromagnet, the dynamical exponent is z = 4. In this case, with the propagator 



the rescaled effective potential reads 



*•<*■"->- yljil+Wt-.- (42) 



V ctf («) = -Ai 2 + isw 4 - (1 + u 2 ) 7 / 2 . (43) 



Higher order terms need to be included at large u to maintain stability. When the <j> field is critical, the third term 
is of order cf) 7 , which is irrelevant. When the <f> field is massive but light, there will also be a term of order (j) 5 which 
is again irrelevant. For general 0, in order for u = to be a local minimum, we need to have A < —7/2. In this case, 
Vl B {u) — has only one positive solution. Thus we have a second-order quantum phase transition. 



Fluctuations in d — 2 and d = 1 



We can calculate the fluctuation-induced effective potential in other dimensions in the same way as above. For 
d = 2,z = 1, and also for d = 1, z = 2, with the rescaled field defined by u 2 = |^-, the rescaled effective potential is 
of the form 

V ctf (u) = -Au 2 + ^Bu 4 - (1 + u 2 ) 3 / 2 . (44) 

When the <fi field is critical, the third term becomes of order — 1-0| 3 , of the Halperin-Lubensky-Ma type, thus the 
quantum phase transition is first-order. Generally when A < — 1.5, AB > — 0.5, B(A + B) > —0.25, u = will be 
a local minimum of the rescaled effective potential V e g, and there are two other local minima at nonzero u. Hence 
there is again a first-order quantum phase transition (see Fig. 5(c)). Otherwise there will be a second-order phase 
transition. 

The effective potential in the case with d = 2, z = 2, and d = l,z = 3 turns out be of the same form as that 
of d = 3, z = 1, as expected from the fact that both cases have the same effective dimension d + z = 4. The case 
d = 2, z = 3 is the same as d — 3, z — 2. 

For d = 1, z = 1, the effective potential takes the form 

V e s(u) = -Au 2 + ^Bu 4 - (1 + u 2 ) ln(l + w 2 ), (45) 

which leads to a first-order phase transition for B < 1 (see Fig. 5(d)). The third term reduces to ip 2 \nip when <fi is 
critical. In this case the quantum phase transition is always first order for any positive value of B. 
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Summary of the classical + quantum cases 

In the table below, we list the most dangerous terms generated from integrating out the fluctuating fields. The 
second row in the table corresponds to the case where (p is critical or massless, and the third row has <fr massive. 



d + z 


2 


3 


4 


5 


6 


7 


massless 


ip 2 Imp 


iP 3 


ip A In ip 




ip e hup 




massive 


(ip 2 + 1) In ip 


ip 3 + ip 


ip 2 Imp 




ip^luip 


V- 5 



One can clearly see that in the massless case, the fluctuations are irrelevant when d + z > 5, while in the massive 
case, they are irrelevant for d + z 7. Otherwise the second-order quantum phase transition can be driven to first 
order. The order of the correction is readily understood from the general structure of the integrals. With effective 
dimension d + z, in the massless case one has 6V/8ip ~ ip J d d+z k(l/k 2 ). Since k 2 ~ ip 2 , this gives the correct power 
SV ~ ip d+z . Replacing gip 2 by gip 2 + and then carrying out the expansion in £ -2 / 1 gip 2 , one gets for the massive 
case a reduction by 2 in the power. We also notice the even/odd effect in the effective potential: for d+z even, 
there are logarithmic corrections. The case d + z = 4 can be easily understood, as the system is in the upper critical 
dimension, and logarithmic corrections are expected. We still do not have a simple intuitive understanding of the 
logarithm for d + z = 2,6. 



TWO FLUCTUATING FIELDS 



We consider in this section the case where the two coupled quantum fields are both fluctuating substantially. The 
partition function now becomes 



Jvip(r,r) f-D<f>{T,T)exp(-S^-S4 



~>4>4>) 



(46) 



We will use RG equations to determine the phase diagram of this system. When there is no stable fixed point, or the 
initial parameters lie outside the basin of attraction of the stable fixed points, the flow trajectories will show runaway 
behavior, which implies a first-order phase transition (72T[76j |. The spin-density- wave transitions in some cuprates and 
pnictides fall in this category 77-9(|. 









f SC 



X 



at^m SSS> ^ 









X 



FIG. 6: Illustration of the fluctuation-induced first-order phase transition in the case of two quantum fields. Here for concreteness 
we consider the antiferromagnetic order and superconducting order. At low temperatures, the phase transitions may become 
first order (the thick vertical lines), due to fluctuations. 



We have considered in the previous sections coupling two single component fields, having in mind that this simplified 
model captures the main physics of competing orders. However, wc will sec below that when the quantum fluctuations 
of both fields are taken into account, the number of components of the order parameters do play important roles. 
So from now on we consider explicitly a ni-component vector field if) and a ^-component vector field coupled 
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-ai\ 



2 + \PiM 4 + l\d^\ 2 



together. When both helds have dynamical exponent z = 1, the action reads 

5^ = J d d rdr 
S^ = J d d rdr 

=9 I d d vdT\^\ 2 \(}>\ 2 



-a 2 |0| 2 + -/3 2 |0| 4 + - 



\d»<t>\ 2 



(47) 



where fj, = 0, 1, • • ■ , <i. This quantum mechanical problem is equivalent to a classical problem in one higher dimension. 
Then one can follow the standard procedure of RG: hrst decompose the action into the fast-moving part, the slow- 
moving part and the coupling between them. The Green's functions are G$ = l/(— 2ai + k 2 + uj 2 ) and = 
l/(— 2a 2 + k 2 + uj 2 ). The relevant vertices are ji\^ 2 ^ 2 pP2^ s ^ 2 pg^ 2 ^ 2 p9^ 2 <^ B ^gi'si>ii>s^f- To simplify the notation 
we rescale the momentum and frequency according to k — > k/A,w — > ui/A, so that they lie in the interval [0, 1]. The 
control parameters and couplings are rescaled according to 0:1,2 — > ai^A 2 , ^x,i Pi,2h?~ d ,g — > gA 3 ~ d . Afterwards 
we integrate out the fast modes with the rescaled momentum and frequency in the range 1]. Finally, we rescale 
the momentum and frequency back to the interval [0, 1], thus k — > &k, ui — > bui, and the fields are rescaled accordingly 
with ijj -> M^/V, </> ~> & (d ~ 1)/2 </>. Using an e-expansion, where e = 3 — d, one obtains the set of RG equations to 
one-loop order, 

^ =2a, - -^[(rii + 2)13,(1 + 2a,) + n jg (l + 2 aj )\, 



8tt 2 ' 

— =< '''-^[{ni + ml 



dl 

dl 
dl 



n 3 g 



(48) 



--9 e 



4tt 2 



[(ni + 2)A + (n a + 2)/3 2 + Ag] 



with index i,j — 1, 2, and i ^= j. These equations are actually more general than considered above. They also apply 
to generic models where two fields with the same dynamical exponent z are coupled together. Generally one has 
e = 4 — d — z, thus a quantum mechanical model with dynamical exponent z is equivalent to a classical model in 
dimension d + z. 




FIG. 7: (Color online) Plot of the RG trajectories in the j3i —@2 plane for two quantum fields with the same dynamical exponent 
below the upper critical dimension. Here we have chosen e = 4 — d — 2 = 0.1. The RG trajectories have been projected onto 
a constant g plane with g = g* , and g* the value of the coupling strength at the stable fixed point, (a) corresponds to the 
case m = U2 = 1, where the fixed point is at /3* = /?J = g* — 4-7r 2 e/(ni + 712 + 8) ~ 0.3948. (b) corresponds to the case 
ni = 2,n 2 = 3, where the fixed point is at (f}*, PZ,9*) = 47r 2 e(0. 0905, 0.0847, 0.0536) ~ (0.3573,0.3344,0.2116). In both cases 
we found that, above some curve (the dashed lines), the RG trajectories flow to the corresponding stable fixed point, while 
below this curve, the RG trajectories show runaway behavior. 

It is known that the above equations have six fixed points [9l|, four of which have the two fields decoupled, i.e., 
g* = 0. They are the Gaussian- Gaussian point at (/3*,/3|) = (0,0), the Heisenberg-Gaussian point at (/S*,^!) = 
(47r 2 e/(rii + 8),0), the Gaussian-Heisenberg point at (y8*,/3|) = (0,47T 2 e/(n2 + 8)), and the decoupled Heisenberg- 
Heisenberg point at 03*,/3f) = (47r 2 e/(ni + 8),47r 2 e/(n 2 + 8)). The isotropic Heisenberg fixed point is at = 
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g* = 47r 2 e/(ni + n 2 + 8), a* = oi* 2 — e(n± + n 2 + 2)/4(ni + n 2 + 8). Finally there is the biconical fixed point with 
generally unequal values of /3*, f3 2 , and g*. In the case, with n\ = n 2 = 1, this is at (/3J", f3 2 ,g*) = 27r 2 e/9(l, 1, 3). For 
Tii = 2,n 2 = 3, one has (Pl,P%,g*) = 47r 2 e(0.0905, 0.0847, 0.0536). 

We find that there is always just one stable fixed point for d + z < 4, below the upper critical dimension [9l|. The 
isotropic Heisenberg fixed point is stable when n\ + n 2 < n c = 4 — 2e + 0(e 2 ), the biconical fixed point is stable when 
n c < rii + n 2 < 16 — nxn 2 /2 + O(e), and when nin 2 + 2(ni + n 2 ) > 32 + 0(e), the decoupled Heisenberg-Heisenberg 
point is the stable one. When the initial parameters are not in the basin of attraction of the stable fixed point, one 
obtains runaway flow, strongly suggestive of a first-order phase transition. Consider for example n% = 2, n 2 = 3, 
where the biconical fixed point is stable. For two critical points not too separated, that is, \a\ — a 2 \ not too large, 
when g > \fp\p2 the RG flow shows runaway behavior, and one gets a first-order quantum phase transition. The 
corresponding classical problem has been discussed in [92[ . We notice the difference from the case with two competing 
classical fields, where one also obtains the same condition for the couplings 7 > V '(3i/3 2 in order to have a first-order 
phase transition. There, the two ordered phases are required to overlap in the absence of the coupling, in other words, 
one needs to have x\ < x 2 . However, in the quantum mechanical case we are considering here, this is not necessary. 
We plot in Fig. 7 the RG trajectories for two cases (a) n\ = n<z = 1 and (b) n\ — 2,113 — 3, where in both cases, 
below some curve, runaway behavior in the RG trajectories is found. 




FIG. 8: (Color online) Plot of the RG trajectories in the fii — P2 plane for two quantum fields with the same dynamical exponent 
in and above the upper critical dimension. The RG trajectories have been projected to a constant g plane. And we have chosen 
ni = 712 = 3. (a) corresponds to the case exactly at the critical dimension with e = 4 — d — z = 0. In this case there is only 
one fixed point with /3* = /3| = g* = 0, the Gaussian fixed point, which is unstable. We found runaway flows everywhere, (b) 
corresponds to the case above the critical dimension, where the Gaussian fixed point is the stable one. Here we have chosen 
e = 4 — d — z = —0.1. We found, below some curve (the dashed line), that the RG trajectories show runaway behavior. 



When d+ z = 4, all the other fixed points coalesce with the Gaussian point, forming an unstable fixed point, thus 
leading to a first-order phase transition (see Fig. 8(a) ). A similar model with an extra coupling and n\ = n 2 = 3 
has been discussed by Qi And Xu [59| . where runaway flows were also identified. Another similar problem with 
d = 2, z — 2 was studied by Millis recently [60|], where a fluctuation- induced first-order quantum phase transition was 
shown to occur. We also notice that in some situations, including fluctuations of the order parameter itself may drive 
the supposed-to-be first-order transitions to second order for both classical and quantum phase transitions 93j-|96[. 

For d+ z > 4, the stabilities are interchanged. The Gaussian fixed point becomes the most stable one. So the basin 
of attraction of the stable fixed point changes. We found numerically that for a given coupling strength g, in the 
Pi — P2 plane, the RG trajectories show runaway behavior when the initial points lie below some curve (see Fig. 8(b)). 
That is, when the coupling between the two fields is strong enough, the QPTs become first order. Just above these 
curves, we found that the RG trajectories will enter the domain with negative /3i or negative f3 2 , and then converge 
to the Gaussian fixed point. For f3i,p 2 large enough, the RG trajectories just converge to the Gaussian fixed point 
without entering the negative domain. 



Competing orders with different dynamical exponents 



We consider next coupling a z — 1 field to another field with dynamical exponent z = z\ 2. To our knowledge, 
such models of two competing order parameters with different dynamical exponents have not been studied previously. 
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The action now takes the form 

S,p = 

S$ = 



d d kduj 
d d rdr 



2 



7i M 



2 fc z i- 2 



d'vdr-p^W 



~* 2 \<t>\ 2 + ~p 2 \cj>\ A 



1 



(49) 



=g / d d rdr|-0| 2 |0| 2 



The new parameter 71 has dimension [71] — L 1 z , and its one-loop RG equation is simply 



(hi i 

— = {z-l) ll . 



(50) 



The Green's function for the ip field becomes 
parameters are modified accordingly, 



l/(-2ai + k 2 + 7i|w|/fc Zl " 2 ). The RG equations for the other 



^1 =2oi - — (m + 2)01 (In 2 + 2a x ) - rWi™2.«(2 + 2a 2 ), 
at 7T7i 

=2a 2 - r> d+1 (n 2 + 2)/3 2 (2 + 2a 2 ) - -^-n l5 (ln2 + 2c*i), 
at 7T7i 



^1 =e/3l _ £^ (ni + 8 )^ _ 2^ +1 n 2ff 2 , 
dl irji 

^ =e/3 2 - 2Q d+1 (n 2 + 8)/3 2 - ^W, 

dl 7T7l 



2fl, 



(51) 



rf 5 ( m d ( . , a on / . x fl fi d 27iln7i+7r 

— =5 e ni + 2)Px - 2^ d+ i(n 2 + 2 )(3 2 ~ 8- — — s .9 

at \ 7T7i 2n 1 + 7i 

where e = 3 — d and 17^ = 27r d / 2 /(27r) d r[d/2] is the volume of the d-dimensional unit sphere. The derivation of the 
above RG equations is included in the appendix. We notice from the above procedure that when the two fields have 
the same dynamical exponent z > 1, one can rescale the couplings to j3\ = /?i/7,/? 2 = p 2 /j,g = 5/7, and these new 
parameters satisfy the RG equations (|48|) with e = 4 — d — z. 

The presence of two different dynamical exponents obviously complicates the problem. It is generally expected 
that the modes with a larger dynamical exponent dominates the specific heat of the system, since they have a large 
phase space, while the modes with a smaller dynamical exponent may produce infrared singularities, since they have a 
smaller upper critical dimension [97] . In the absence of the coupling between the two fields, we have the RG equations 



^L ={ A-d-z)h-—{ni+S)Pl 
al it 



(52) 



dfa 
dl 



--(3-d)p 2 -2n d+1 (n 2 +8)p 2 



For d = 3, P2 is marginal with an unstable fixed point, while pi is irrelevant and its Gaussian fixed point is stable. 

Generally, for z > 1, if the initial value of 71 is nonzero, the absolute value of 71 will increase exponentially. The 
RG equation for /3 2 becomes independent of other parameters, 



^ = e p 2 -2Q d+1 {n 2 + 8)p 2 . 

We are interested in the case e = 0, for which /3 2 is readily solved to be 

1 



A(0 



P2 1 + 2Q 4 (n 2 + 8) (/-Z cr )' 



(53) 



(54) 



with p 2 taken at the crossover scale l cr at which the P\ term begins to dominate the g 2 term. Only the sign of p 2 
matters. If (3 2 > 0, as I increases, /3 2 will decay to zero, flowing to its Gaussian fixed point. From the simplified RG 
equations for g, 

dg 
dl 



-2n 4 (n 2 + 2) 9 p 2 , 



(55) 
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fa Pi 




FIG. 9: (Color online) Plot of the RG trajectories in the /3i — @2 plane for two coupled quantum fields with different dynamical 
exponents. The RG trajectories have been projected to a constant g plane with g = 1. We have chosen the spatial dimension 
to be d = 3, the dynamical exponents Z\ = 2, Zi = 1 and the number of field components n\ = ri2 = 3. (a) shows the RG 
trajectories originating from the region below the dashed line, which flow to negative /3i or negative fa regions, (b) shows the 
RG trajectories originating from the region above the dashed line, and those flow to the stable points on the positive axes of 
fa, the location of which is sensitive to the initial value of the parameters.. 

one can see that with a lower power in fa, g drops to zero even more quickly than fa. Taking fa as quasi-static when 
considering the evolution of g, one notices that g decays exponentially as g(l) ~ exp(— 2^4(712 + 2)fal). So dfa/dl 
also decays exponentially, and before fa goes to zero, fa already stabilizes to a finite value fa, which depends on the 
initial value of fa. Actually from the simplified RG equations for fa,fa,g with I/71 set to zero, one can see directly 
that the fixed points are at fa = g* — 0, with fa any real number: we have a line of fixed points. When fa > 0, there 
will be a second-order phase transition. When fa < 0, the transition becomes first order (see Fig. 9). 

If fa < 0, the absolute value of fa will increase without bound. Subsequently g and fa also diverge, leading to 
runaway flows. 

CONCLUSION 

Quantum criticality in the presence of competing interactions is an important guiding concept that allows us to 
organize a framework for emergent states near QCPs. Here we investigated the stability of a quantum critical point 
in the presence of competing orders. We focused on a simple quadratic-quadratic interaction, where coupling between 
two competing phases is assumed to be of g"4> 2 4> 2 form. We find that QCPs are often unstable and transform into 
first order lines of transitions. The detailed scenario on how the instability develops depends on the precise nature of 
the competing interactions, dynamical exponents and strength of the coupling. The general trend we observe is that 
competing interactions, be they classical or quantum, often lead to the instability of QCPs. This instability in fact 
always occurs, in the cases we have investigated, if the coupling g is strong enough. We thus conclude that breakdown 
of QCPs is a ubiquitous phenomenon. The magnitude of the specific heat jump in some first order transitions (the 
classical + classical case) is of the same order as the specific heat released in a second order transition and these 
first order transitions are strong, and not weakly first order as found in Halpcrin-Lubensky-Ma. An immediate 
consequence of this breakdown is that we can expect spatially modulated inhomogeneous phases to be present near 
QCPs, given their propensity to turn into first order transitions. The wide likelihood identified here of first order 
transitions preempting a QCP leads us to anticipate the nucleation and metastability phenomena associated with 
such transitions [981 ] . Additionally, proximity to first-order transitions makes auxiliary fields (e. g. magnetic field, 
strain) and disorder very important over substantial parameter regions [99j . 

The broad similarities we pointed out between QCPs and AdS/CFT models offers an interesting possibility that in 
fact AdS models arc also spatially inhomogeneous. More detailed analysis that allows breakdown of scaling, specific 
for AdS/CFT is suggested. 

We derived the renormalization group equations for two coupled order parameters with different dynamical expo- 
nents. We found that there are a line of fixed points, which is quite different from the case where two o rder parameters 
have the same dynamical exponent. Very recently, there have appeared some interesting reports (97L 1 1 00j ] investigat- 
ing the eff ects of the presence of two order parameters with different dynamical exponents near the Pomeranchuk 
instability |l0l| . as examples of multiscale quantum criticality. It would be interesting to see how the presence of 
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two different dynamical exponents, and the coupling between the corresponding order parameters, affect the scaling 
of resist ivity , especially whether a linear- resistivity is possible, overcoming the "no-go" theorem for single parameter 



scaling [102 1 



In this paper, we have confined ourselves to the framework of Hertz-Millis-Moriya 68-7(|, considering only the 
interplay of bosonic order parameters. It would also be interesting to study the electronic instabilities, to see whether 
the superconducting instabilities and Pomeranchuk instabilities are enhanced in fermionic quantum critical states. 
Fermi liquids, even wi th r epulsive interactions, are unstable towards forming a superconducting state, due to the 
Kohn-Luttinger effect 103j | resulting from the presence of a sharp Fermi surface. Fo r the fermionic quantum critical 
states, the momentum distribution function may have only higher order singularities [104j. It would be interesting to 
check whether the Kohn-Luttinger effect is still active in this case. 
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APPENDIX: DERIVATION OF THE RG EQUATIONS FOR TWO FIELDS WITH DIFFERENT 

DYNAMICAL EXPONENTS 



In this appendix, w e wil l derive the RG equations of two competing orders with different dynamical exponents. We 
follow the notation of |l05j . Our starting point is the action ([49]) . First we count the dimensions of the field operators 
and all the parameters: 



M - [r] =L, 
[k] = [u] =L-\ 

M = [4>] =l^ 2 , 

[ai] = [as] =L~ 2 , 
[Pi] = m =L d -\ 
[g] =L d -\ 

[7i] 



(56) 



Then we decompose the action into slow and fast modes. The action for the slow modes reads 



e( s ) fiw / , fc2 _l T 1 M 



Si s) = / d d vdr 



- a2 |0 s | 2 + i/3 2 |^| 4 + i|d^| 2 



4 S i =3 / d d rdT\xp s \ 2 \4> s 



(57) 
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Since we will only consider RG to one-loop order, the interaction terms in the fast modes, the contraction of which 
leads to second-order diagrams, can be ignored. Thus we obtain the action for the fast modes, 



c(/) _<?(/) + c(/) 



s 



(/) 



d kdu —at 



d d rdr 



-a 2 \4> f \ 



from which one can easily identify the Green's functions as 



7i M 
2 fc 2 i- : 



(58) 



g{M=- 



-2ai + A: 2 + 7i^ 
Si i 



(59) 



-2a 2 + k 2 + u) 2 

The coupling between the slow modes and fast modes takes the form 



d d vdr 



(3/3x^^^ + 3^44^^) +9\^s\ 2 \^f\ 2 + g^f\ 2 \ct>s\ 2 + Ws-ipf)^s ■ </>/) 

ijkl 



(60) 

with the tensor F ijk i = \{8 i:j 8 k i + Si k 3 3 i + SuSjk). 

Now we can integrate out the fast modes and see how the different parameters change accordingly. The effective 
action of the slow modes is determined by 



exp 



-5'. 



exp 



exp 



2\ con 



-{s c ) f + -(Si)} 



(61) 



In the S 2 term we take a connected average, thus the superscript "con". The coefficients in the RG equations 
will depend on the different renormalization schemes. Here we will u se the procedure that is most convenient for 
the problem at hand, similar in spirit to what was outlined in [l06| . We integrate over the momentum interval 
A/6 < k < A, which after rescaling k —> k/A, gives b^ 1 < k < 1. The frequency part is more complicated. We 
will introduce a cutoff when it is necessary, otherwise just integrate over the whole range — oo < u> < oo. The main 
reason for us to choose this RG scheme is that in calculating the third correction to the coupling g, the two internal 
lines come from order parameters with different dynamical exponents, thus the two frequencies scale differently with 
momentum, and this RG scheme offers a simple and self-consistent treatment of the cutoffs. 
71 does not receive corrections up to first-order. 



First Order Corrections to at\ 




FIG. 10: One-loop diagrams contributing to the first order correction of ay,. The solid lines represent the ip fields, and the 
dashed lines represent the <j> fields. The external lines are slow modes, and the internal lines are fast modes. 
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Two terms in the action (|60| contribute to the first-order corrections of a\. The coupling tp 2 ip 2 leads to the 



correction 



SVS[a 1 ]=3fk y £F ij H ^|^(4(q0^(-q0>/^r^(q)^(-q). (62) 



ijkl 

Using the identity, 



one obtains 



y^-Pjifei = ni t^ &kh (63) 



6^S[a 1 ] = (n 1 +2)^ [ du f dfcfc^ 1 — -i / J^^, a (q) . ^.(-q). (64) 

2tt 7 A-i -2ai + fc 2 + 7i J s (27r) a+1 

Assuming that the V 1 field is near its critical point, thus a\ is a small parameter, the Green's function can be expanded 
in terms of — 1tt\. The correction term can be written as 

6<»S[<*i] = (»i + 2 )^^( J i + 2 «i J 2) J (^pTT^(q) ■ tf.(-q), (65) 
where we have defined the series of functions 

r 1 1 
duj / dkk d ~ x - - Tr . (66) 

J ^ (* 2 +7i^) 

Let us first calculate I\. The frequency integral requires a cutoff. From dimensional analysis, we choose to integrate 
over the region — 1 < 71 u> < 1, and obtain the result 

h = ~f In (i±^V (67) 



7i ib-i V fc2:i 

To proceed further, we are required to specify the dimension and dynamical exponent. Consider d = 3, z\ =2, where 
one has 



37i 



1 + b~ 2 

In 2 - &~ 3 ln ( ) + 2(1 - 6" 1 ) - 2 arctan 1 + 2 arctanfc -1 



(68) 



Expanded to first order in (1 — b 1 ), it is simply 

21n2 



For d = 2, Z\ = 2, we obtain 



/ 1 = ^rz(i_ 6 -i). ( 69 ) 
71 



Ji = — [2 In 2 - (1 + b- 2 ) ln(l + b- 2 ) + In 6~ 2 1 , (70) 
7i 



which leads to the same result when expanded to first order in (1 — b 1 ). This result can also be obtained more 
crudely by setting fc = 1 in the integrand of (|67|) . I2 can be calculated similarly, with the result 

I 2 = -(l-b- 1 ). (71) 
7i 

So the one-loop correction to a\ coming from the coupling ip 2 ip 2 is 

S^S[ax} = (m + 2)^-^-(l - 6- 1 )(ln2 + 2a x ) f -^-^(q) . ^ s (-q), (72) 
7T7i J s (2ir) d+L 
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We next calculate contributions from the coupling tp^cfyj, which takes the form 

5^S[ ai }=g^ jkl jf ^^(^(q')^(-q')>^^^^(q)^(-q), (73) 

with F[- kl = SijSki- So we have simply the identity 

Y, F kki =n 2 S k i, (74) 

i 

which gives 

s m s[ai] . „ M g jtof dkk^ — + \ 2 - -, / ^*.(q) - ».(-q). (75) 



Defining the new set of functions 



''^hC^WTW* (76) 



one obtains 

S^S[ ai ] = n 2 g^(I[ + 2a 2 I>) jf ^^.(q) ■ ^.(-q). (77) 
Here we integrate over frequencies in the range — oo < w < oo, and get 

I[ = 7T / dfcfc d ~ 2 , (78) 

which is, to first order in (1 — fo -1 ), 

I'^Al-b- 1 ). (79) 

Similarly for 1 2 we have 



j£ = 1 f dkk <i-\ 
l Jb- 1 



thus 



(80) 



I' 2 = l{l-b~ l ). (81) 



Near d — 3, one has fid/4 — fid+i- Grouping all these together, we obtain the second term in the correction to ct\ as 

5^S[a{\ = n 2 gn d+1 (l - 6" 1 )(2 + 2a 2 ) J A^L^q) ■ tf.(-q). (82) 

First Order Corrections to «2 

The calculation of the first order corrections to a 2 is quite similar to that of ot\. There are again two terms 
contributing. The coupling V'/^s gives rise to a term of the form 



/Jd+l-y f f ]d+l n 

j^WMW-J)) J s ^p^(q)^(-q), (83) 



t ijkl I 

ijkl 

Summing over the field indices, 

J2 F uM=ni\4>s\ 2 , (84) 
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FIG. 11: One-loop diagrams contributing to the first order correction of 0,2- The solid lines represent the fields, and the 
dashed lines represent the <j> fields. The external lines are slow modes, and the internal lines are fast modes. 



we obtain 



which can be expanded as 



The result is 



ijkl 

We first sum over the field indices, 



£Wta = ^l^| 2 , (89) 



3 

resulting in 



Expanding to first order in a 2 , one has 



and the final result is 



^JL0 8 (q) . a (-q) 



First Order Corrections to /3i 



(86) 



S«S[a 2 ] - m ff -^-(l - b^ln 2 + 2ai) / -f-3_0 s (q) • s (-q). (87) 

7T7l J s (27Tj a+1 

The other term comes from the coupling <t> 2 f<j) 2 s - It has the form 

^-^<4(q')^(-q')) j( («,)#(-<,). (88) 



). (90) 



5P)S[a 2 ] = (n 2 + 2)/^^ +2a 2 / 2 )^^^0 s (q) • s (-q), (91) 



(92) 



The first order correction to Pi comes from two one- loop diagrams, one with two internal "0/ lines, the other with 
two <pf lines. The dependence of the internal lines on the external momenta and frequencies can be ignored here, 
since they arc of higher order. 
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FIG. 12: One-loop diagrams contributing to the first order correction of /3i. The solid lines represent the V fields, and the 
dashed lines represent the <j> fields. The external lines are slow modes, and the internal lines are fast modes. 

The first term with ipf internal lines is of the form 

—-^ (V} 1 (q>) 2 (-q')> W (q')V? (-q')> (g3) 
x J d d rdTii) k s 1 ip k s 2 ip l s 1 il) l s 2 + 2 permutations. 

Using the identity, 

E F ijkihFijk 2 i 2 = g[(ni+ 4:)S klh S k2 i 2 +24 1 fc 2 (5i 1 ( 2 +25^5^], (94) 

combined with the 2 other permutations of the external lines, the part containing the field component indices can be 
simplified as 

E E^^^V^^ 2 + 2 permutations = ^±^ s | 4 . (95) 

kik 2 hl 2 ij 

Thus the first correction to (5\ reads 
which is, to leading order of ai, 

= -(m + 8)^^/2 y d d rdr|t/> s | 4 . (97) 

Substituting the explicit expression for I2, we get the result 

SWSfa] = -(m + 8)ft 2 — (1 - / d d rrfr|^ s | 4 . (98) 
I'll 7 

The second term has two <pf internal lines, and takes the form 

zii 2 hh k lk2 hi 2 JiV*) K^) 

x J d d rdT^ k 8 1 \p k s 2 \p l s 1 ^ l 8 2 + 2 permutations. 



The part with the field component indices gives 



E E ^mA^V^VW + 2 permutations - n 2 |t/> s | 4 , (100) 

kik 2 hl 2 ij 
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which further leads to the result 

WSVk] = J toj^ ^ d - \_ 2a2+ \ 2+oj2)2 J (101) 

To leading order in a 2 , it is 

5^ Si/3,] = -n29 2 ^I' 2 J d d vdr\^ s \\ (102) 



or more explicitly, 



«J< 2 >S[ft] = -n 2 g 2 n d+1 (l - b- 1 ) J d^dr]^. (103) 



First Order Corrections to 




FIG. 13: One-loop diagrams contributing to the first order correction of /3 2 . The solid lines represent the ip fields, and the 
dashed lines represent the cj> fields. The external lines are slow modes, and the internal lines are fast modes. 

The first order correction to the j3 2 term also comes from two diagrams. The first one has two ipf internal lines, 
and is of the form 

svs\fh] = -g > £ £ F> ijikih F! M /^^{^(^(-^{^(^(-q')) , m 

»i»2jij2 fciMi^ f 

x J d d vdT(j) k 8 1 <l) k s 2 (j) l s 1 <l)^ + 2 permutations. 
Summing over different field components, where one has 

E J^ F Urh F ljk 2 lJs 1 <P': 2 <P l s<P l s 2 + 2 permutations = m|0 s | 4 , (105) 

k 1 k 2 hh ij 

the first correction to the f3 2 term is 

S (1) S[/3 2 } = -n l9 2 — [ duj [ dkk d - Y n / d d rdT\<t> s \ 4 . (106) 

To first order in a i, it is simply 

SWSfa] = ~ni9 2 ^l2 J d d vdr\<^ s \\ (107) 



which can be written as 



5 W S [p 2 } = -n 1 g 2 ^-(l-b- 1 ) f d d rdT\cf> s \ 4 . 



(108) 
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The second diagram contains two (f>f internal lines, thus the correction reads 

<5 (2) S[/? 2 ]=-(3/3 2 ) 2 J2 E F iMF i2j2k2l2 /'^^(^(q')^(-q'))(^ 1 (q')^ 2 (-q')) (1Q9) 

k 1 k 2 hh iiiihh S y > V I 

x J d d TdT^ 1 <t>g 2 4> l ^(t> l s 2 + 2 permutations. 
The summation over the field indices gives 

E Y^ F V^ F iih*l*<fc<fc&& + 2 permutations = ^^|^| 4 . (110) 

k\k 2 l\l 2 ij 

Thus the second correction to /3 2 reads 

wsm = -(n 2+8) ^| / ^f h _^ { _ 2a2+ \ 2+uj2)2 J d^m\ (in) 

When the </> field is near its critical point, the above expression can be simplified to be 

8 {2) S[p 2 ] = -(n 2 + 8)/3 2 2 ^/ 2 | d<W|0. 



which is 



First Order Corrections to g 





(112) 



<5 (2) S[/y = -(n 2 + 8)/3|f7 d+1 (l - 6- 1 ) y d<W|</> s | 4 . (113) 



FIG. 14: One-loop diagrams contributing to the first order correction of g. The solid lines represent the ip fields, and the dashed 
lines represent the <f> fields. The external lines are slow modes, and the internal lines are fast modes. 

There are three diagrams contributing to the first order corrections of the coupling g between the squares of the 
two fields. The first diagram has two ipf fields as internal lines. This term takes the form 

s {1) s[ 9 ] = - - x 2 x 2(3/3^5 E E F nh^ F Uk 2 i 2 I T^r^CqO^C-q'))^^^)^^-^)) 

k 1 k 2 l 1 l 2 iii2jij2 f (114) 



f, i ,h jk 2 



x / d d vdr^^c 



The 1/2 comes from (1/2)5^, and the two 2 factors come from the expansion in S% and the number of contractions 
in (ipfipf(x)ipfil>f(y)). We first sum over the field indices, 

E E^^^b^^ 1 ^ 2 ^^^!^! 2 !^! 2 , (us) 

kik 2 hh ij 



and then substitute the Green's functions, 



= -2(3 l9 (ni + 2)^ / dcj f dkk^ 1 — -± r— f d<W|V> s | 2 |0, ' 2 

2lr J A- 1 (-2a 1 + fc 2 + 7l J^) 2 J 

Keeping only the leading order term, 

6^S\g] = -2p 1 g(n 1 + 2)^I 2 J d d rdr|t/> s | 2 |</> s | 2 , 



one arrives at the result, 



5^S[g] = -2p ig ( ni + 2)^(1 - b- 1 ) f d d vd T \^ s \ 2 \<}> s 

7T71 J 



5^S[g]=--x2x2(3l3 2 )g E £ F 4l jlklh F> 2 jM I —-^ (q')# (-q')><<# (q')<^ 2 (-q' 

fcifco/i/o itini-i in f 



The internal lines of the second diagram are two <f>f fields. The corresponding correction term is now 

d d+i„> 

t , i 1 j 1 k 1 l 1 -t , i 2 j 2 k 2 l 2 I 

k 1 k 2 hh *i»2jij2 

x J d d vdr^^ k s^ l s 2 - 
The summation over field indices is similar to the first term, 

E E^m^^I 1 ^ = ^i*.i a i*.i a - 

kik 2 hh ij 

Thus the correction to the action is also similar, 

8^S[g] = -2p 2 g(n 2 + 2)^- J dw J dkk^ 1 ^ + 1 fc2 + ^ 2)2 | d d rdr|^| 2 |0 s | 2 , 
which is, to leading order in a 2 , 

5^S[g] = -2p 2 g(n 2 + 2)^I' 2 J d d rdr|t/> s | 2 |0 s | 2 , 



or 



S^S[g] - -2f3 2 g(n 2 + 2)fl d+1 (l - b' 1 ) J d d vd T \^ 8 \ 2 \cj> s 



The third diagram has one 4>f internal line, and one ipf internal line. The correction takes the form 

d d+i / 



8^S[g]=-\{Agf £ E F 'i M K M J |^(^(q')^(-q'))(^(q')^/(-q')> 

ktk 2 l t l 2 i\ii3\3i f 



With the summation 



we obtain 



k 1 k 2 hh ij 

6^S[g} = -8g 2 ^ [ du t dkk^ 1 — L= ? ^ / <fW|t/> s | 2 |0 s | 2 . 

LJ 2^7 y 6 -! -2a 2 + k 2 +w 2 -2 ai + k 2 +^ lT ^L J l^ll^i 



Assuming both fields are near their critical points, the above equation is approximately 



S^S[g] = -8g 2 ^l" J d d rdrW s \ 2 \<j>. 



S I 7 
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with the new function I" defined as 



1 



1 



I" = du dkk^ 1 - 



(128) 



As mentioned before, here in our RG scheme, frequency is integrated over the whole real axes. In the two propagators, 
frequency scales differently with momentum. For the ip field, jiu> ~ k Zl ; for the (j) field, u> ~ k. And a finite cut-off in 
frequency would lead to inconsistencies for such cases with miscellaneous dynamical exponents. Here the frequency 
integral gives 



dkk 



d+zi— 3 



1 



b- L 



?k 2 + k 2z i 



To leading order in (1 — b ), we have 



I" 



11 



1 



221—2 



-71 In ■ 



k 2z ^ 



nk 



21-1 



1+7? 

This leads to the third term in the correction to the g term 



(l-6- 1 )(2 7 i ln 7 i+7r). 



2i l+7f 



(129) 



(130) 



(131) 



Rescaling of the Parameters 

Now combining all the above results for the corrections of the different parameters, and carrying out the rescaling 



we obtain the RG equations 



r — s-r/fo, 
t ->r/6, 
^ _ >6 (rf-i)/a^ j 



(132) 



7i l 7x, 



-ai ->6 2 
-a 2 -^6 2 



2 
2 



n 



ai + (m + 2)A — (1 - 6- x )(ln2 + 2a x ) + n 2 <?Oi+i(l - JT 1 )^ + 2a 2 ) 



7T7l 



-a 2 + (n 2 + 2)(3 2 Q d+1 (l - b- 1 )(2 + 2a 2 ) + n ig —(l - 6~ 1 )(ln2 + 2a x ) 

7T71 



^ - (m + 8)^-^(1 - 6- 1 ) - n 2ff 2 n d+ i(l - 6- 1 ) 

2 7T7l 

h _ m ff 2 -^(l - - (n 2 + 8)/3|^ + i(l - b- 1 ) 
g - 2Mm + 2)-^(l - 6- 1 ) - 2$ 2 g(n 2 + 2)fi d+ i(l - b^) 8g 2 ^ + * (1 - 

7T7l Z7T 1 + 7f 



(133) 



the differential form of which has been presented in equations (j5TJl l5Tj) . 
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